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A new powerful method to test the stability of the replica symmetric spin glass phase is 
proposed by introducing a replicon generator function g(v). Exact symmetry arguments are 
used to prove that its extremum is proportional to the inverse spin glass susceptibility. By 
the idea of independent droplet excitations a scaling form for g(v) can be derived, whereas 
it can be exactly computed in the mean field Sherrington-Kirkpatrick model. It is shown by 
a first order perturbative treatment that the replica symmetric phase is unstable down to 
dimensions d < 6, and the mean field scaling function proves to be very robust. Although 
replica symmetry breaking is escalating for decreasing dimensionality, a mechanism caused 
by the infrared divergent replicon propagator may destroy the mean field picture at some 
low enough dimension. 
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I. INTRODUCTION 

The theory of spin glasses has by now become more than three decades old, 1 even if one consid- 
eres the archetypal model by Edwards and Anderson [6] as its beginning. Nowadays spin glasses 
are regarded as the prototypes of complex systems showing highly nontrivial equilibrium properties 
(massless glassy phase without long-range order, nonergodic Gibbs state with a complicated free 
energy surface, lack of self-averaging of some quantities, etc. . . . ) on the one hand, and a very 
slow aging dynamics with diverging relaxation times on the other hand. The state of affairs in the 
understanding of basic issues of the equilibrium glassy phase is still controversial, as two rather 
different theories have emerged for the low-temperature phase of the short-range Ising spin glass: 
One of these two schools extends the idea of replica symmetry breaking (RSB), invented by Parisi 



(see 



and references therein) for the mean field version of the model — i.e. for the Sherrington- 



Khkpat rick (SK) model [jj — , to the short-range case in the physical dimensions. I I • oilier one. 



'Electronic address: temtam@helios.elte.hu 

1 For reviews at different stages of the story, see P, 0, S 01 an< i the collection of papers in [f| about different aspects 
of spin glass theory. 



the so called scaling or droplet picture [8|, |9|, LLC]) [ill] , claims — at least in its most radical version 



- that RSB is restricted to the SK model which, with its infinite coordination number, can be 
relevant only in infinite dimension (d = oo), or for long-range models in finite d. Before displaying 
a list of the most important conflicting statements, we present the model and give some definitions. 

The Hamiltonian of the Edwards- Anderson (EA) model on a d-dimensional hypercubic lattice 
with the newest neighbour q nenehed interactions J i/s taken .ndependently tan a 8ymm et„c 
Gaussian distribution is [6] 

(ij) » 



with N Ising spins sitting on the lattice sites, and N — > oo in the thermodynamic limit. (The 
second term contains the external magnetic field H, although it is zero throughout this paper.) The 
distribution Pj{q) of the site overlaps qss' = j? SiSj, while picking out the spin configurations 
S and S' independently by the Gibbs-measures with the same Hj's, is usually considered as the 
order parameter function. The maximum value of the Pj(g)'s support is the EA order parameter 
gEAi although originally qj = X^(£i) 2 = I dqqPj(q) was proposed as qea p, which generally 
depends on the realization J. Here and in subsequent definitions (...) means thermal average by 
the Gibbs measure ~ e —Hj/kT ^ w h ereas 777 stands for averaging over the quenched disorder. We 



have P(q) = Pj(q). The following three 2-site (zero momentum) correlation functions are useful 
to define: 



N 



G2 = ^E^->^>^->-^ 2 



(i) 



G 3 = ^E^) 2 ^) 2 - Ar ^ 2 - 



N 



Two linear combinations of these are especially important to characterize the spin glass phase {§]: 
G R = Gt - 2G 2 + G 3 = ^ ((WW = Xea (2) 

ij 

and 

G L = Gi-4G 2 + 3G 3 . 

The subscripts R and L are for "replicon" and "longitudinal", and Gr = xea is the famous 
Edwards- Anderson (or spin glass) susceptibility whose inverse is the so called replicon mass 2m\ = 
XeA' which is used throughout this paper to test the stability of the spin glass phase. 
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In the following list, it is attempted to review the most important conflicting statements provided 
by the two rival theories: 

1. The Gibbs measure is decomposed into an infinite number of pure states, not related by any 
symmetry, in RSB. The droplet picture is like the ferromagnet: there is one pair of pure 
states related by spin inversion. As a result, a trivial probability distribution of the overlaps 
between equilibrium spin configurations arises, with a pair of delta functions at i(?EA- A 
nontrivial anc 
Parisi theory 

delta function for self-overlap, cjea is, however, self-averaging in both theories. 



non-self-averaging [i.e. Pj(q) ^ P(q)] overlap distribution follows from the 
2I, with a continuous part for a range of overlaps smaller than (/ea plus the 



2. The ergodic components, or pure states, in the Parisi theory are ultrametrically organized 

a. 

3. The average system-size low-energy excitations, when modelling in a finite size system, are 
proportional to N®/ d , > being the stiffness exponent, and the surface of the compact 
cluster of flipped spins, i.e. of the "droplet", is a fractal with d s < d in the droplet picture 
Pillll- On the other hand, the same type of exitations are thought to be of order unity, and 
they are space filling (d s = d) in RSB. 

4. The external magnetic field H destroys the spin glass transition in the droplet picture, 
whereas the spin glass phase is bounded by the Almeida-Thouless (AT) line Hp^(T) in RSB 
[jjj]. Along the AT line xea of ([2]) is infinite, that is the replicon mass is zero and starts to 
become negative, signalling the end of the replica symmetric (RS) paramagnetic phase, and 
the onset of RSB. 

5. There is one feature both theories share, namely the spin glass susceptibility xea is infinite, 
in the thermodynamic limit, in the whole glassy phase, i.e. it is a marginally stable massless 
phase. The details, however, already differ, as the divergence of xea with N is 



Xea ~ < 



N for RSB, 

(3) 

N i-e/d droplet picture, see [8, 9j. 



The discussion of the above points may be started by considering "replica equivalence," a set 
of identities between joint distributions of the site (or standard) overlaps, and rigorously proved 
for the SK model |l j - Q, Its first derivation was based on the ultrametric property of the 



Parisi solution of the SK model [lfl, it is, however, more generic than ultrametricity [17]. Applying 
replica equivalence to the susceptibilities in ([I]), we can easily derive for the leading O(N) term 
[see also ©]: 



Go 



5*. 



and 



XEA 



3 3 



dqq 2 P(q) 



dqqP(q) 



This formula is in complete accordance with point 5 above for both scenarios. It is quite interesting, 
and not completely understood, that the ratio 1 : ^ : g for the leading terms of the G's, which are 
0{N 1 - Q l d ), in the droplet case remains true, although not because of replica equivalence. 

Replica equivalence has recently been rigorously proven for a large class of models, including 
the EA model, for the link overlap q l ss , = ^ ^iSj S^Sj 



study 



18 



19, 



20l | . A very recent numerical 



21[ for the 3-dimensional model provided strong support that the two overlaps (site and 



link) are equivalent in the description of the spin glass phase of the EA model. Considering also 
the fact that replica equivalence is a very plausible property when deriving it in the replicated 
system (with the replica number n going to zero) [171 ]. it can now be believed that the identities of 
replica equivalence are also true for the standard site overlap in the finite dimensional Ising spin 
glass. As both scenarios, RSB and droplet, fulfill replica equivalence, it can not, however, be used 
to distinguish them. 2 

To review the huge amount of numerical work accumulated in the last decades is impossible, the 
reader is referred to Refs. 24], |25|] for extensive lists of papers. Numerical simulations work on finite 



systems and try to extrapolate to the infinite N limit. Due to the extremely long relaxation times, 
notwithstanding the huge evolution of computer capacities, only relatively small systems have been 
simulated. That resulted in a fluctuating support for the two scenarios. Ref. [24] considered almost 
all topics of the above list, finding evidences for replica equivalence and RSB. A very recent paper 
261 ] for the three-dimensional ±J model in zero external field found excellent agreement with 
ultrametricity — the lattice size was relatively large: 20 3 — , which is a unique characteristic of 
RSB. The issue of the low-energy excitations (point 3) is rather controversial: Refs. 



27, 



2J found 



system size O(l) energy excitations, but with d s < d suggesting that the link overlap distribution 
is trivial, whereas the standard one is not (TNT picture). This scenario, however, was ruled out 



by 



211 ] . Zero temperature methods were also used to test the survival of the spin glass state in 



2 It has emerged as an exact statement in the literature [§2, [23] that the distribution of the site overlaps for any 
finite-dimensional EA model must be self- averaging — P (q) = Pj{q) — in case of the infinite system. Putting 
together this thing with replica equivalence, one can deduce a trivial overlap distribution, i.e. a droplet scenario 
in any d. Marinari et al in Ref. [241 ] argued that the Pj{q) defined by Newman and Stein is different from that 
appearing in the replica equivalence identities. Nevertheless, one might have the feeling that this issue has not 
been completely settled by now. 
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a nonzero external field 



29(; the somewhat positive answer obtained is in favour of RSB. On the 



n 



other hand, the completely different treatment in [30[] concluded the lack of a transition in field; 
though the magnetic field was not homogenous in that investigation, but chosen from a random 
Gaussian distribution. Finally point 5 is rather difficult to check numerically as the exponents in 
([3]) are very close to each other, due to the smallness of Q/d. 

The replica method introduced by Edwards and Anderson [6], is a powerful, though admittedly 
unphysical, tool to handle the quenched state by eliminating the inhomogenities caused by ran- 
domness. The effective (replicated) Hamiltonian thus obtained can be represented by a field theory 
using common procedures known from the study of ordinary systems. The advantage of field theory 
comes from its flexibility of application in any space dimension d and for a variable range of the 
interaction. Standard methods, like systematic perturbation expansions and the renormalization 
group (RG), are also available. In the first calculation by replica field theory of Harris et al 31] 



the fixed point value and the critical exponents were obtained in leading order in e = 6 — d. In two 
subsequent papers 



32 



331 ] the ordered phase below T c , and especially the stability of the RS state 
in finite d was studied: the mean field, AT-like, instability was found to survive perturbatively. Af- 
ter Parisi's solution of the SK model appeared and proved to be (marginally) stable 34|, one-loop 



studies based on the ultrametric RSB ansatz were published (see 35], and references therein), the 



severe infrared problems were handled in the framework of a scaling theory 



361 ] . In the meantime, 



a search for finding the (presumably zero-temperature) fixed-point of a finite-dimensional AT line 



took place 



33, 



38]. 



In field theory, detection of a glassy phase consistent with the droplet picture would mean 
finding a marginally stable RS state in zero external field. The developments enumerated in the 



previous paragraph were against droplet theory in finite dimensions. In Ref. [321 ] a procedure with 
temporarily breaking the replica symmetry by a special — non-Parisi — type of RSB (the "two- 
packet" RSB), and eventually restoring RS led to a massless RS phase. We will comment this 
result in Section IVIIl The issue has been revived in two recent papers \3i 



vve 



Alike the solution of the fully connected SK model, field theory provides results in the ther- 

attice EA model is, however, 



411 ] . i.e. averaging out short- 



modynamic limit, i.e. for N — » oo. A direct comparison with the 
impossible due to the use of the concept of the minimum length 
range fluctuations. As a result, the precise dependence of the parameters of the field theoretical 
Lagrangian, the bare couplings, on the control parameters (temperature, external magnetic field, 
etc.) of the original model is lost. This fact may be extremely important when trying to distin- 



guish an AT line from a droplet-like phase 



42 ]: both have the generic replica symmetry with a 



6 



ass susceptibility (a zero replicon mass). 



421 ] . of using exact symmetry arguments 



nonzero EA order parameter, and even an infinite spin gl, 
For that reason, we follow the ideas, introduced in Ref. 
to parametrize the bare couplings in zero external field. Our arguments are based on the only 
assumption that the Legendre transformed free energy is an invariant of the extra symmetry of 
the high-temperature phase, see Eq. ([7]) below. An important result of the present paper is the 
summation of the infinite terms which were argued in [42| to contribute to the replicon mass 2mi: 
there is a complete cancellation of these terms leading to a negative mass, i.e. to an unstable RS 
phase for 6 < d < 8 and d < 6. 

The outline of the paper is as follows: In Section [TT] the generic replica symmetric model is 
introduced, and the extra symmetry of the replicated paramagnet of the high-temperature phase 
presented. Our basic object, the replicon generator g(v), is defined, and its properties are derived 
by exact symmetry considerations. The full g(v) is computed for the SK case in Section IIIII The 
picture of independent droplets is used to derive the scaling form of higher order susceptibilities 
in Section IIVI A somewhat heuristic argument makes it possible to produce a scaling form of 
g(v ) for the droplet scenario when the replica number n is nonzero. Exact Ward identities, which 
are necessary to construct the free propagators, are displayed in Section [V] whereas the one-loop 
calculation of g{v) is presented in Section IVIl The qualitatively different regimes — i.e. d > 8, 
6 < d < 8 and d < 6 — are discussed in separate subsections. Section IVIII is for a summary of the 
results, and it also contains our conclusions. Notations, definitions and computational details are 
left to the series of appendices, amongst them the block diagonalization of a generic matrix with 
the "two-packet" RSB is in Appendix [El 



II. THE SYMMETRY OF THE SPIN GLASS FIELD THEORY IN ZERO MAGNETIC 

FIELD 

The Ising spin glass on a hypercubic <i-dimensional lattice — with the quenched exchange 
interactions taken independently from a Gaussian distribution with zero mean — can be represented 
(see [43y and references therein) by a replica field theory where the classical fields <p a ^ have the 
Artfa replica indices a, » = 1 . . . n, and ^ = with z e r „ d.ngonai f-. In the spirit of the 
replica trick [6j, the special limit n — > of this n(n — l)/2 component field theory provides the 
physical results. Unlike in the Parisi theory, the spin glass phase is almost ergodic in the droplet 
scenario: the two minimal free-energy phases are related by the global inversion of the spins - 
just like in a ferromagnet. In a replicated picture this gives rise to an RS theory with the generic 
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Lagrangian (displayed here in momentum space): 



c = \ E f {\p 2 + O E <K +^E <^- 7 P +^E ^p^-p 

a/3 a/37 a/37<5 



2 

p 



+ £ r , (4) 



with the interaction term containing all the higher order invariants I- together with the single 
first order one, i.e. 

c = -Ni I - E E' E (s) 

a/3 fc=3 iV 2 Pi i 

The number of spins N goes to infinity in the thermodynamic limit, its explicit indication ensures 

- -(k) 

that the bare couplings, h and the g)- s, are independent of N. Momentum conservation is 

(k) 

understood in the primed sum over the momentums pi . . . pk- The I- s are invariant under the 
global transformation <f>"^f = 4>p aPl3 , where P is any permutation of the n replicas. The three 
invariants of the kinetic term {k = 2) and that of the linear one are explicitly displayed in @ 

and ([5]), while the 8 cubic and 23 quartic ones are summarized in Appendix [A] 3 Besides using 

Jk) ■ ■ ■ ~ - 

-I ii / 1 1 1 / ii' i / ■ /• Ti-ii t i'wi ininii mil' lu'Ari i 'i / 1 1 -w inriL' I, is- + lo^v order couplings h ?Ti's 



Qj for a generic k, we retain the commonly used notations for t 



u;'s (for cubic) and u's (for quartic). As it was pointed out in 42 ]. an ambiguity in assigning the 



bare coupling constants to a given physical state occurs due to the freedom to offset the zero- 
momentum fields, <^p£ — * ^p=o — l eavm g an the irreducible vertices unaltered while the 
bare couplings changing. By this way, we can always tune the one-point function to zero, hence 
avoiding the appearence of any "tadpole" insertions in a perturbative treatment. We adopt this 
definition of the bare coupling constants in what follows. (In this case, h is generally nonzero even 
if the "physical" magnetic field is switched off.) 

The generic replica symmetric theory presented in the previous paragraph, however, is not con- 
fined to the description of the droplet-like spin glass phase: the nonglassy phase in an external 
magnetic field (even above the transition temperature) is also represented by the generic RS La- 
grangian of Eqs. (|4l5p . Still, as a spin glass is always an infinitely correlated state, the criterium 
of a massless mode (the so called replicon one) seemingly may distinguish it from the magnetized 
paramagnet. At that point, however, we must realize that the boundary of a possible RSB phase 
fits the same criterium: along the AT line we have replica symmetric massless thermodynamic 
states. 



3 The barred notation for the bare coupling constans is used throughout the paper to distinguish them from their 
exact counterparts. The choice of the set of coupling constants of the generic replica symmetric Lagrangian is 
obviously not unique: couplings with the lowercase notation defined above are the coefficients in front of the 
invariants. To make the presentation as clear as possible, it will be necessary later to introduce another set with 
the capital letter notation (mi vs. Mi, for instance); the linear relationships between these two sets are displayed 
— for the quadratic, cubic and quartic vertices — in Appendix iBl 
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The underlying problem we face when trying to distinguish an AT line from the droplet-like 
spin glass phase is the indirect relationship between the bare couplings of the field theory and the 
physical control parameters (temperature and magnetic field, for instance). To find the thermal 
route of the system, i.e. the dependence of the bare couplings on the reduced temperature Tc ^ T 
while keeping the external magnetic field zero, we turn to exact symmetry arguments. The "true" 
paramagnet, i.e. the high temperature phase in zero magnetic field, has the extra symmetry that 
the Edwards- Anderson order parameter q is zero. (The magnetic field makes q nonzero even in the 
nonglassy phase.) As a result, those couplings in (JH) and ([5]) which belong to invariants with at 
least one replica occuring an odd number of times (like m,2, rn^ and h) are exactly zero. In addition 

n 

to the permutational invariance, this results in the extra symmetry published in [421 ] . and for the 
sake of self-containness, it is explained here again: The n replicas are divided arbitrarily into two 
groups with p and n — p replicas, p being a free parameter between and n. The Lagrangian is 
invariant under the global transformation 



b^f for a and 3 in the same group, 

(6) 

-0p^ for a and 3 in different groups. 



This transformation is obviously orthogonal in the n(n — l)/2-dimensional space of the field com- 
ponents, with the diagonal transformation matrix O a p^s = (— l)° n ^ +1 g. We introduced here 
the overlap a n 3 defined asan/? = 0(aH/3 = l)ifa and 3 are in different groups (if a and 3 
are in the same group). 

The symmetry of the Lagrangian is reflected in the invariance of the Legendre-transformed free 
energy F: 

F(q' af3 ) = F{q ap ). (7) 

The stationary condition -§^—^ = provides the symmetrical solution q' a ^ = q a p = for the 
paramagnet, while in the droplet-like phase the extra symmetry of © is lost leading to q a n = q ^ 
represented by the generic RS theory in (|4|) and ([5]) . In this way, we can understand the paramagnet 
to spin glass transition symmetry breaking one. 

In [421 ] the invariance property of ([7]) was combined with the freedom to choose the continuous 
parameter p of the transformation to derive a set of Ward-like identities between exact irreducible 
vertex functions (derivatives of F evaluated at stationarity). Here we focus on the stability of 
the low-temperature RS phase in zero magnetic field, and for that reason we define the generator 
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function g(v) of replicon-like vertices as follows: 

oo 

9(v)=Y,Ti Gik+1)vk > (8) 

where 4 



fe=2 



G {k+1) = E' A; = 1,2,...; and by definition G« = G$. (9) 

(a 1 /3l) i ...,( Qfe /3 fc ) 

The prime above is to indicate the restriction ctj n = for all z = 1, . . . , k in the summations, 
and we choose a n /? = too, hence G'(' c+1 ) no more depends on a and /3 in an RS state (it does, 
however, on n and p). The generic definition of the exact vertices is 



, (fc) _ d fc (F/AQ 

. . . dq akl 3 k 



G axPx,...,u k P k ~ »„ a_ ( 10 ) 



evaluated at an RS stationary point. For k = 2, 3 and 4 we insist — for traditional reasons - 
to use the notations —M, W and U, respectively, similarly to their lowercase counterparts, see 
the note below Eqs. (|4)5p . The vertices defined in (jlOH take only a limited number of different 
values, expressing in this way the RS symmetry; in fact they may be numbered in accord with 
the numbering of the set with the lowercase notation. All of the relevant definitions are left to 
Appendix (Aj while the linear relationships between the two sets are displayed in Appendix |B] for 
the quadratic (k = 2), cubic (k = 3) and quartic (k = 4) cases. (See also the footnote [3j) 

Differentiating both sides of Eq. (|7|), and exploiting the orthogonality and diagonality of the 
transformation in ([6]), the following rule is easily derived for the transformation of the vertices: 

G ai/3i,...,a k /3 k = ^aiPi,aiPi ■ ■ ■ O ak p k ^ ak j3 k G^fc^^fa, (H) 

inferring immediately 

G (*+i)' = (_l)*+i G (k+i). (12) 



As it was explained in 42|, the transformation q a p — * q' a/3 becomes infinitesimal for q, n, p <C 1, 
hence the left-hand side of (fl~2j) can be Taylor-expanded around the RS state q a p = q. From 
and (fTUj) . and using q' a g — q a p = unless a fl (3 = 0, when it is — 2q, we can easily derive: 

oo 1 

G (k+l)' = J2 i G( k+1+ V(-2q) 1 . (13) 
z=o 



4 Throughout the paper, we use different notations for summations over distinct pairs of replicas, — i an< l 

(a/3) a</3 

unrestricted sums, = Yl X)i as m ^ ne invariants in Eq. Q. 

q/3 a 
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Comparing ([12]) and (JTSJ, 



oo 1 



/! 

z=i 

is concluded, valid for = 0, 1, ... . It is useful to analyze the following cases separately: 

• For k = 0, using (151), -2 GO-) = G( 2 ) x (-2*7) + g(—2q) is obtained. Stationarity renders 
the left-hand side zero. Here we give some details of the derivation of G^ in terms of the 
masses, intending in this way to facilitate the reader to check more complicated formulae for 
G^ and G^ displayed in Appendix O 

= ^ M a ^ s = M 1 + [(p-l) + (n-p-l)}M 2 + (p-l)(n-p-l)M 3 ; (14) 

( 7 <5) 

for the notations used here see Appendix [A] This formula and the similar ones displayed in 
Appendix O are much simpler when expressed in terms of the lower case vertices, though 
the derivation uses evidently the upper case set. By Eq. (|Bip of Appendix [Bl we arrive to 
the important result 

2m\ + nm,2 + 4p(n — p)ms = ^ evaluated at v = —2q. (15) 

• For k = 1 the right-hand side is just g'(—2q), i.e. 

g '(v) =0 at v = -2q. (16) 

• The generic case, k > 2, is easily understood by observing that the right-hand side is 
gW(-2q) - G^ k+l \ while G^ k+l ^ = g( k \0) providing 

{ _ 1) k + l g (k) {v = Q)=g (k) {v = _ 2q y 

The above conditions are serious restrictions for the RS field theory representing a spin glass 
phase in zero magnetic field. For a given q (which is used here — instead of the reduced temperature 
- to parametrize the spin glass phase below T c ), the generator g(v) must have an extremum at 
v = —2q, Eq. (|16|) . while the higher derivatives satisfy (|17|) . As 2m\ is the dangerous (so called 
"replicon") mass responsible for the instability of the RS state [12|], stability requires g(v = —2q) < 
in the limit n — > 0. (Equality means marginal stability, i.e. a massless spin glass phase, just what 
droplet theory suggests [a, Q].) These results were derived by generic symmetry arguments, hence 
they must be equally valid for systems described by mean field theory and for short-range spin 
glasses in physical dimensions. All these cases are reviewed in subsequent sections. 
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III. CALCULATION OF THE REPLICON GENERATOR IN MEAN FIELD THEORY 



An explicit expression for F(q a p) is available in mean field theory, working it out either on 
the fully connected lattice (Sherrington-Kirkpatrick model or as the long-range limit on the 



(i-dimensional hypercubic lattice 



4^: 



jv^ = ^E<4- ln Jr } exp ( £ q^FSP) . 



(18) 



(a/3) v (a/3) 

As we are in zero magnetic field, the only control parameter is Q ~ (fcT) -2 . It is now useful to 
introduce the modified Legendre-transformed free energy F by 



F(v a/3 ) = F(q + v a p), 



(19) 



q being the RS stationary point of F. Obviously F is stationary for v a p = 0, and all the vertices 
in Eq. ()10p can be computed equally well substituting F by F and the q's by v's. (The invariance 
property of (|7|) is, however, lost for F.) After Taylor-expanding F and making the restriction 
v afi = v f° r a H P = 0, otherwise it is zero, we can express g(v) by the now one-variable function 
F(v): 



9(v) 



d 



{F/N). 



p(n — p) dv 

The second term above is the harder to compute. As it is common in mean field replica calculations, 
the Gaussian integral representation method can eliminate the spin traces in (|18j) providing for 
p(n — p) — > 0: 



giy) = 

f V y f V x f V x i tanh(- v /g + v y + y/—v x) tanh(- v /g + v y + \J — v x') cosh n (yjq + v y + \J—v x) 

J v y fv x j v x , cosh n (^ny y + x) 



tanh 2 (^/gy) - 1 + (n - 2) tanh 2 (^/gy) - (n - 1) tanh 4 (^/g?/) v. (20) 



The shorthand notation 



tanh fe (^y) 



J X> v tanh fc (v^y) codify) 



with 



e 2 



V 



J^cosh n (^y) J a J V2n 

was used here. The explicit temperature dependence has been eliminated by the stationary condi- 



tion G l q = tanh 2 (y / gy), and g(v) is thus parametrized entirely by the order parameter q. 
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We are now in the position to check the generic properties of g(v) deduced in the preceeding 
section. It is relatively easy to see that g(v) is quadratic for small v. Although it is somewhat 
harder, a straightforward application of the following two lemmas 5 can prove Eqs. (|16p and (11T|> : 



1. 



d 
dv 



Jv y J V x j V x , h{z)f 2 (z') = jv y jv x j V x , fi(z)f^z') 



and 

2. 

jv y jv x j V x , h{z)f 2 (z') = j T> z f 1 ( y /qz)f 2 (- y /qz) for v = -2q. 

fi and /2 above are two arbitrary functions, while z = \/q + v y+^—v x and z 1 = y/q + v y+^—v x'. 

Near T c , g(v) gets into the simple scaling form g(v) = q s g(v/q); the scaling function results 
from expanding the integrands of (|20p for q <C 1 and v/q = 0(1): 

g{x) = -^^x 2 {3 + x), p(n-p)^0. (21) 

The following remarks are appropriate here: 

• g{— 2) > 0, showing by (fT5|) the instability of the RS mean field spin glass phase just below 
T c . 

• The leading, or scaling, term of g(y) — and only that — satisfies the following condition: 

g(v = -2q) = ^q 2 g"(v = 0). (22) 

We shall see later that this property of the replicon generator is very robust, as it persists for 
the short-range system below the upper critical dimension 6. 6 On the other hand, a simple 
one term Ward-identity follows from (|22p for the replicon mass, as one can see from f)15[) and 
422): 

2m 1 = ^ ( _ = 2 ^ ] " = ~\ q G^ = -\q W2 , P ,n^0. (23) 



5 The first one can be easily proven by a sequence of partial integration, while a change of variables to z, z' and 
y' = y/yj—v in the second lemma and setting v = — 2q renders the integral over y' proportional to 5(z + z'), leading 
us immediately to the end of the proof. 

6 To be more precise, the one-loop calculations in Section IVTl will show that Eq. (I22|l is true for any n when d > 8, 
whereas for 6 < d < 8 and d 5; 6 it is restricted to the case n = 0. 
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[The expression for in terms of the w 7 s is displayed in Appendix [Cl see Eq. ()Cip ,] This 
formula is valid whenever the leading contribution to g(v) for n — ► and q — * satisfies 
(122p . As it connects the replicon mass {2m\) and the replicon-like cubic vertex (1U2), it is 
obviously important in testing stability. We shall therefore return to it later. 

• The leading term of g{v) in Eq. (|21h can be deduced directly from the truncation of the 
infinite series ([8]) as 

g(v) = 1 -G^v 2 + 1 -G^v 3 . (24) 

This can be seen by a direct evaluation of G^ and G^ from the definitions Q, ()10p and 
using the mean field free energy (jl8[) . This truncated form of the scaling part of the replicon 
generator is, however, more general, as it remains valid for the finite-dimensional system 
when d > 8 (see later). Starting from (|24|) . Eq. (|22|) trivially follows, together with all of its 
consequences. 

From the generic properties of g derived in the preceding section and from the mean field analysis 
above, one may get the feeling that the replicon mass [which is proportional to g(v = —2q)] is either 
negative or positive, but it is difficult to imagine how marginal stability may emerge. This, however, 
belongs to the essence of the droplet picture, and in the next subsection we try to construct g(v) 
— or rather its scaling term — using the ideas of the theory of free droplets. 



IV. THE REPLICON GENERATOR IN THE DROPLET PICTURE 



The droplet picture is a low-temperature scaling theory developed for the original d-dimensional 

1,3, Im 



lattice system 



111 ] . As a basic inference of this theory, the spin glass state is now restricted 



to zero external magnetic field, and it is in fact a critical surface attracted by a zero-temperature 
fixed point, immediately suggesting that temperature is now a (dangerously) irrelevant parameter. 
The whole phase below T c is massless, as at least one correlation length, namely that of the replicon 
mode, is infinite. Somewhat surprisingly, there is a single independent exponent belonging to the 
zero-temperature fixed point — the stiffness exponent 9 — , and scaling of relevant and irrelevant 
variables, and also of several observables (correlation functions, for instance), can be described by 
means of it. 

Translating the results of the droplet theory to the language of replicated field theory is straight- 
forward by using the Gaussian integral representation, called Hubbard-Stratonovich transforma- 
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tion, and then eliminating short ranged fluctuations by the introduction of effective coupling con- 
stants and truncation of the momentum-dependent mass term (see Ref. [43]); a Lagrangian like 
that in Eqs. (jlj) and ([5]) is obtained. As a short cut to translate different quantities, we propose that 
whenever a pair of replicated Ising spins SfS^ occurs (i being a lattice site), it should be replaced 
by the field </>°^ ; in this way, Ising averages are turned to averages with the measure e~ c . As a 
pedagogical example, let us consider the momentum-dependent Edwards-Anderson susceptibility 
defined by 

xea(p) = ^E gip(r ^ ri) (Mhw, 

y 

where (...) and 777 mean thermal average with the Ising measure and quenched average over the 
disorder, respectively. Applying the replica trick and introducing the fields, we get 

— {kT) b («^ p > - 2(^0%) + «Vl* p » = (kT) b cW(p); 

(p) being the replicon component of the 2-point connected correlation function in the replicated 
field theory. (To avoid complicated notations, we stick to (...) for thermal averages with different 
measures: in the first line it is taken with the effective replicated Ising Hamiltonian resulting 
from the replica trick, whereas a field-theoretic average with the weight e _-C is understood in the 
second one.) The emergence of the factor (kT) b above results from the transformation. Similar 
temperature-dependent factors with some simple exponents always appear whenever spin averages 
are represented by field averages. Near T c they are absolutely irrelevant, we must, however, keep 
them at low temperature to ensure the correct temperature scaling. (Nevertheless, we shall see 
that the specific values of b and also a, which is introduced later, will not influence our basic 
conclusion.) We know from Refs. 0, [ijj that Xea(p) ~ (kT)p~ d+e yielding 

C (2) (p) ~ (kTf- b p- d+6 '. (25) 

We can apply the scaling formula 

around the zero-temperature fixed point. \i\ is a relevant scaling field vanishing in the spin glass 
phase (e.g. the external magnetic field) with Ai > exponent, while \i2 ~ {kT) a is the dangerously 
irrelevant one with A2 < 0. The scaling index of the replicon correlation function is —2 + t/r, 
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and the reader is reminded that, in principle, there are two different 77's (t/r and 1%) at a generic 
RS fixed point when n = 0. On the one hand, we can now exploit the fact that must be 
independent of /ii when it goes to zero, and the correct momentum-dependence must be ensured. 
On the other hand, droplet theory suggests A2 = —a6. From these, we can conclude: 

C^{p)^(kT) 2 - zm ^p- d+e . (26) 

Comparing ([25]) and ([26]) . we get 

VR = 2-d + b6. (27) 

The replicon mass [2mi; see the table for the masses in Appendix [A] and (|B1|) ] is — by Dyson's 



equation 44] — just the zero momentum limit of (p) _1 , and therefore obviously zero for any 
T in the spin glass phase. 

The momentum-dependence of the higher order correlation functions C^(p), k > 2, can be 



deduced by applying the single droplet exitation picture of Fisher and Huse 
Edwards-Anderson susceptibilities. As an example, for k = 3 we have: 



9, 



111 ] for generalized 



i.^ e *(pir l+ p 2 r J+ p 3 r fc ) ((C^.^) _ (^(S^) - (S J )(S l S k ) - (S^SiSj) + 2{Si){Sj){S k ))* 
ijk 

— > 5 Pl+P2+P AkTfc^(p), b' = \b, 

where momentum conservation is explicitly shown, and only a single p scale is considered. C^(p) 
has the structure of 2w2, see (|Blh and the table for the cubic vertices in Appendix |Aj and it is fully 
replicon-like [43]. The lengthy formula for C^(p) is omitted here; the reader can reproduce it by 
replacing a cubic vertex W a( g )7 j j/U/ in 2w2 by the average y/~N (0pf <pp S 2 (j>ps) , 7 and also symmetrizing 
with respect to the momentums. In the spirit of the droplet theory, we can conclude the momentum- 
dependence of the cubic and even higher order replicon correlators as 

C {k) (p) ~p e -( k - 1 ) d : k = 2,3,... . 

It follows from simple scaling arguments that the replicon correlators, which are by definition 
independent of N (see footnote [7|) , have the scaling index d+ |(?7r — 2 — d), and for fi2 <C fi\ X2 ^ Xl 
we have 

C^(p) Ml [*H-*foR-2-<i)]/2Ai (^/^V*!)"* cW^iA^ (28) 



The \/N factor here is to make C'^ 3 \p) become independent of N, i.e. of order unity. For the same reason, we 
define C (fc) (p) with the prefactor AT' -1 . 
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where the one-variable scaling function C^ k \x) ~ x 6—(k—l)d £ or x _^ ^ whereas — as the replicon 
correlation length is finite for fj,% ^ — it goes to a constant in the opposite limit, is nonzero, 
expressing the dangerously invariant nature of fi2, and it can be computed by taking the limit 
fii — > 0, while H2 and p finite, of C^ k \p) in ([28]) : 

Kk X 2 = k(^f-l + fj -0=(±b-l)0, (29) 

where the last equality follows from (|27p . We shall need the scaling of the zero-momentum corre- 
lation functions when approaching the spin glass state for a given temperature. It can be easily 
deduced from (1281) and (1291): 

C( fc )(p = 0)~ ^-(k-i)d]/x^ kT>Q fixed (3Q) 

It is pointed out in Appendix O that G (fc) , Eqs. ^ and CE]), is a fully replicon vertex in the 
spin glass limit, and therefore — applying the usual connection between vertices and connected 
correlation functions [44| — it scales like C <yk \p = 0) x C^ 2 \p = 0)~ k . Using ([30]) . this gives 

G i.k) ^ ^-{k-m/x^ kT>0 fixed (31) 

This equation, which is valid for n = 0, shows that g(v) of the spin glass in the droplet theory 
is highly singular for [i\ — > 0, i.e. for zero magnetic field. This is in full contrast with mean field 
theory, see Eq. (|2"Tj) . Nevertheless, (JUT]) can be used to guess the n-dependence of g(v) in the droplet 
picture in the following way: Although Eq. (|31|) was deduced by renormalization group arguments, 
it merely expresses the assumption that an infinitezimal magnetic field drastically destroys the 
droplet picture by introducing a finite correlation length. Regarding the common observations 

n 

(see [43] for instance) that the role of the replica number is alike, i.e. it makes the replicon mode 
more massive, we can assume tentatively the following form for the n-dependence of g{v) in zero 
magnetic field: 

g (v) = n dp/e g(v/n p ). (32) 

This equation is obtained by the formal replacement of [i\ by n in (|3il . and then substituting G (fc) 
into the definition of g(v) in Eq. (|8j). A new independent exponent p was introduced here, and 
- although not explicitly indicated — the function g depends on T too. We shall return to this 
formula later, extending it to the critical scaling region around T c . 
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V. WARD IDENTITIES AND THE FREE PROPAGATOR FOR T < T c 



Some of the identities between exact vertices presented in Ref. [42] follow directly from the 
properties of g(v). Eqs. ([8]) and (fl5|) provide: 

/ 2 \ 

2mi + nm 2 + 4p(n - p)m 3 = - [G {3) - -G (4) q + ...jq, (33) 



whereas it follows from both (116]) and (jT7j) : 

G( 3 ) = g + . . . . (34) 
One more expression arises from the transformation property of a mixed mass component 
M a/3, 7 «5 = - M «/3, 7 5, a H (3 = I and 7 n5 = 
[see (fTT|) ]. summing over 7 n £ = and expanding the left-hand side. We get for p — > 

m 2 = -[(wi + ~w 3 ) +n(-u; 5 -l-we)] <? + ••• • (35) 



This equation has also been displayed in Ref. 



42|. Everywhere in the above formulae, the dots 



have the obvious meaning of higher order vertices multiplied by the appropriate power of q. 

Eqs. (|33p . (|34p and (|35p are now used to construct the bare propagators for a one- loop calculation 
of the leading scaling term of g[v). The following two regimes are separately treated: 

• In the perturbative regime (d > 6) the system is defined by the set of bare couplings com- 
patible with the symmetry of the high-temperature phase: w\; ui, U2, U3, U4; ... . The 
zero- loop limit of (|33|) and (|34p . using also ([C2]) of Appendix O gives 

2 1 
2mi + nm 2 = --n 2 g 2 , 2m 3 = --(«i + 2u 4 ) q 2 ; 

whereas ([35]) provides 

m 2 = —w\ q- 

Neglecting the q 2 term in the scaling limit, the free propagators assume the diagonalized 



form 



43] 



, with rR = n(wiq), r^ = 2(w\q) and n, = (2 — n)(wiq). (36) 



p2 _|_ r 



[The multiplicities are: n(n — 3)/2 for R (replicon), n — 1 for A (anomalous) and 1 for L 



(longitudinal); for the details, see 
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For d < 6 we enter the nonperturbative regime. Nevertheless, for 6 — d = e <C 1 the fixed 



point wl is small, and the perturbative renormalization group results of Ref. 4a] are usable. 
By the scaling analysis in [42] we have shown that q k ~ 2 ~ wf 2 (w\q)~& , with w^ 2 = O(e) 
and j//3 = 1 + 0(e). The right-hand side of ([33]) and ([35]) are of order e, with the only 
exception of — w% q in Eq. (j35[) . This results in the bare masses 

2fh\ + nrri2 = 0[e{w\q)}, fh^ = 0[e (w\q)\ and f?i2 = — {w\q) [1 + 0(e)}. 



After replacing vb\ by w*, the free propagators remain the same as in (|36p . 
In conclusion, in the one-loop calculation of the next section we can use as bare masses: 

-w\q for d > 6 

2fhi + nm,2 = 0, 7713 = and rri2 

-w\q for d < 6. 



(37) 



VI. THE ONE-LOOP CALCULATION OF 

Our basic formula for g(v) can be derived by a Legendre-transform method; leaving the technical 
details to Appendix [Dl we only quote the result here: 



g(v) = g(v) +H + vY^ (M a p, j5 

(r/S) 



9{v) 



^ v 2 ^2 WaPrfS^v + ^ ^ Uap,ys,nu,pu + • • • , and all (3 = 0. 

[Alike in ([U]), the primed summation means the constraint that the replicas in a summing pair 
belong to different groups.] The second and third terms ensure that g(0) = g'(0) = 0, as it must 
be by the definition in ([8]), whereas g(v) is the replicon generator built up from the bare couplings. 
Due to the choice v a p = v for a n = 0, and otherwise zero, the partition function Z of the 
system with the couplings in (|D4[) depends now on the single variable v. Obviously InZ has the 
"two-packet" RSB structure, nevertheless it serves as the generator of the replicon-like vertices of 
the RS system. Although Eq. (|38|) will be used in a first order perturbative calculation in this 
section, it is quite generic, and may be useful even in a nonperturbative treatment. 
The one-loop contribution to the last term of (|38p is 



1 



d 



p(n — p) dv 



1 

N 



InZ 



with 



(38) 



9i(v) 



1 



1 



2p(n -p) N 



EE 



1 



p 2 + Xj dv 



(39) 
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where Xj is an eigenvalue of M, the bare mass of C displayed in (|D4h . The restriction v a /3 = v 
when a n (3 = 0, while otherwise it is zero, means that M is a matrix with the "two-packet" RSB 
32l . l39| . whose block diagonalization is worked out in Appendix [E] for the most generic case. In 



the scaling limit {q — » and v ~ q), M reduces to the first two terms in (|D4j) . and only W\ is kept, 
as it is the only cubic coupling which is finite at criticality. With the notations of Appendix [El we 
then have 



M 



(i) 



M, 



(i) 



M. 



M 



Mi 



M 



(2) 



M 



'(2) 



M 2 , 



,>(3) _ ^/(3) 



(40) 



M 



M 2 - Wi 



M. 



(i) 



M. 



(2) 



M 



'(2) 



M 



(3) 



M 



'(3) 



M. 



(4) 



Af s . 



'3 — — ^3 — ^"3 — ^"3 — ±V± Z ~ ±vl 3- 

As it turns out from Appendix [El the two replicon classes are independent of v in this approxi- 
mation, and they do not contribute to g%. To simplify the following algebra, and exploiting the 
arbitrariness of p, the symmetric choice p = n — p = n/2 m computing the longitudinal (L) and 
anomalous (A) eigenvalues will be applied. After substituting the mass elements from Eqs. (|40p 
and ([37]) into the L and A blocks of Appendix [El we get 



2w\q 



1 +n(n - 2)(l + 



n 



q — > and v ~ q. (41) 



wiq 



n + 2 



1 ti\ 2 

l + -n(n-4)(l + - 
4 a 



Putting these expressions into (|39p . with the multiplicities taken from Appendix [El the following 
result is obtained for the one-loop scaling contribution: 



9l(v) 



n 



n 



n 



( p 2 +A + )(p 2 + A - ) 4 (p 2 + A + )(p 2 + Aa) 



(42) 



(For ci < 6, w\ must be replaced by w\.) It is obvious from (|4ip and (|42|) that gi satisfies 

- 2q) = -gi(v). 



By a direct calculation of the second and third terms of (|38p . it was checked that g\{v) in 
ensures that the generic property g(v = 0) = g'(v = 0) = is valid up to one-loop order. 

By straightforward, though somewhat lengthy, manipulation of Eqs. (|38p and (|42p — whilst 
using the definitions of the ^-independent masses in (]36p and the f-dependent ones in (I4ip — , the 
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following convenient form for the one-loop leading scaling term of g(v) can be concluded: 

g( v ) = g( v ) + 4(n - 2) W! (w ig f {v/qf (3 + v/q) - ^ 

p 

(n-4) 2 /16 



+ 



-v. 



(p2 + A+)(p 2 + A~)(p 2 + r L )(p2 + rR ) 

+ 8?i(n-2)u)i(u)ig) 5 (?;/g) 2 (2 + t;/g) 
(n - 4) 3 /64 



(p 2 + A+)(p 2 + AX)(p 2 + r A )(p 2 + r R ) J 
(n - 2) 2 

(p 2 + \£)(P 2 + Al)(;p 2 + r L ) 2 (p 2 + r R ) 2 ' (p2 + A + )(;p 2 + A - )(:p 2 + rA) 2 (;p 2 + rR) 



+ 



(43) 



This formula was derived for p = n/2 in the scaling limit g — > and v ^ q. In what follows, Eq. 
(|43p is used to track the behaviour of g(t>), and along with it the question of stability of the RS spin 
glass phase, starting from high dimension (d > 8) and ending below the upper critical dimension 
(d < 6), where the nontrivial fixed point is small, and one-loop results can be translated to the 
leading behaviour of g(v) in an e = 6 — d expansion. 



A. The mean field like high dimensional regime: d > 8 

In this case, the bare replicon generator g(v) as well possesses the properties of Eq. (|17p as the 
exact one. From this and the definition in (|38p . the leading term for q — > and q ~ v follows as 



1 1 

g{v) = - [2u 2 + nu 3 + -n 2 (^i + 2u 4 )] q 3 (v/q) 2 (3 + v/q), (44) 

where (|C2p of AppendixOhas been used for p = n/2, and only the four quartic couplings remaining 
mite at T c are kept. [To reproduce the mean field results of the Sherrington-Kirkpatrick model 



, we must substitute the values u\ = 3, U2 = 2, u 3 = —6 and U4 = into (jM|). We must be 
careful, however, if a direct comparison with Eq. (|2ip is to be made, as it is valid in the limit 
p — > 0, whereas p = n/2 was chosen in this section. Obviously, disregarding the n 2 term in (|44p . 
which comes from p(n — p), leads to complete agreement.] Putting together this result and the 
high momentum leading term of (|43|) . g(v) assumes the scaling form q 3 g(v/q), with 

g(x) = j i [2u 2 + nu 3 + X -n 2 (u x + 2n 4 )] + ~ (n - 2)(n 2 + 8n - 16) w)f ^ ^ -if x + x ) 
I p P J 

A comparison with (|2ip reveals the surprising fact that the scaling function g[x) ~ x 2 (3 + x) 
remains very robust with respect to mean field theory (d = 00), apart from a normalization factor, 
at least up to one-loop order. 
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The replicon mass for finite n can now be easily derived from (|15p and ()35|) : 

2m 1 ^nwiq- (45) 
2g 

For n > 0, there is a temperature domain below T c where the RS phase is stable (2m i > 0). The 
condition 2mi = signals the end of the stable RS phase, defining an Almeida-Thouless line in 
the T-n 46], or equivalently in the q-n plane (instead of the original AT line in the temperature- 



magnetic field plane): 

QAT — C<2 n, (46) 

where the dimension-dependent slope has a finite limit for d — » oo (the zero-loop mean field limit), 
whereas it goes monotonically to zero while d — * 8; see Fig. [TJ A first order perturbative correction 
leads to a shrinking RS phase (or growing RSB phase) when decreasing the space dimension. 

B. The transitional regime: 6 < d < 8 

Although critical behaviour is still governed by the Gaussian fixed point, and standard per- 
turbation expansion remains valid (no accumulating infrared divergences), scaling of g(v) is now 
modified to 

g(v) =w 1 (w 1 q) 2 ~ t/2 g(v/q), e = 6 - d. 



The scaling function g follows from (|43p after sending the cutoff to infinity and rescaling the 
momentum p 2 — ► p 2 /(wiq). (The bare contribution g(v) is subleading for d < 8, and it can be 
neglected in what follows.) The following properties are easily concluded: 

1. For n = 0, the mean field form g ~ x 2 (3 + x) is again recovered, together with its immediate 
aftermaths in Eqs. (|22p and ([23]). The truncated form of (|24p then follows as well. This is, 
however, not at all trivial when d < 8 — as it was pointed out in 42| — , since all the terms 



in ([8]) are now the same order, and a cancellation of the terms higher than cubic order in v 
must occur. This solves the problem raised in Ref. 42J: the infinite sum of its Eq. (11) is 



exactly zero, leading to instability. This conclusion is now reached by the evaluation of the 
g(v) function, and the previous statements are valid, of course, in this first order perturbative 
framework. 

2. A nonanalytic n-dependence develops for d < 8, which is proportional to n 1-£ / 2 . Although 
this nonanalycity is innocuous with respect to the zero n limit, the truncated form (I24p is 
no longer valid for finite n. 
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The boundary of the RS state follows again from (|45p and (|43p . the zero n limit of the propaga- 
tors taken from ([36]) and (I4ip . The AT line in the q-n plane is now nonanalytic, and the tendency 
of a shrinking RS phase persists when approaching the upper critical dimension 6 (see Fig. [T]): 

1 



n 



p 4 (p 2 + 2) 2 



-2<e = 6-d<0. 



(47) 



C. Below the upper critical dimension: d < 6 



,*2 



The nontrivial fixed point [31f] w\ 
dimensions. A scaling analysis in Ref. 
Eq. ([8]) in the scaling limit 



-e/(n — 2) + . . . governs critical behaviour below six 
42j] showed that G^ k+1 ^ q^ 1 ~ u)f (w 1 g) 7//3 , yielding by 



g{v)=w\ (wtqy/^giv/q), 



d<6. 



(48) 



With this definition above, g has a finite limit for e — > 0, and it obviously de pen ds on d and n. 
Renormalization group calculations provide 7//? + 1 = 2 + e/2 + . . . , see 
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451 ]. Our result is 



consistent with this form at the leading order in e, the scaling function is the same as that of 
the transitional regime in the previous subsection, except that the momentum integrals must be 
performed at d = 6. The first property of the previous subsection, concerning the n = behaviour 
of the scaling function, can be entirely taken over, whereas the n-dependence of g, which is harmless 
around d = 6, will be discussed in the conclusions. 

Both terms in Eq. (|45p are now proportional to [w\q)'~ 1 / ^ , and we get, in the limit n 
means of (|43p the replicon mass 

d 6 p 1 



0, by 



n 



+ Q{e 1 



e = 6 - d > and n > 0. (49) 



(2vr) 6 p 4 (p 2 + 2) 5 

We can see from this formula that the stable RS spin glass phase for n small but finite — which 
does exist in mean field theory and even for d > 6 — completely disappears below six dimensions, 
continuing in this way the trend of a shrinking RS phase (see Fig. [T]). 



VII. SUMMARY AND CONCLUSIONS 



The Ising spin glass in zero magnetic field has been studied in this paper, using the field 
theoretic representation of the replicated model and concentrating on the question of stability of 
the replica symmetric glassy phase. Any preconception concerning the RS spin glass phase has 
been avoided, except that the Legendre-transformed free energy preserves the symmetry of the 
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replicated paramagnet, Eqs. ([6]) and (|7|), just as in the case of common phase transitions with 
spontanous symmetry breaking. As a convinient way to study stability, the replicon generator 
function g(v) was introduced. On the one hand, it was exactly computed for the Sherrington- 
Kirkpatrick model, whereas droplet theory ideas were applied to guess its behaviour in case of 
a droplet picture scenario. In this later case, a highly singular behaviour in the spin glass limit 
(n — > 0) was found. We followed the fate of the mean field (d = oo) scenario by a first order 
perturbative calculation down below the upper critical dimension 6, where the fixed point is small 
and the perturbative renormalization group is valid. The following conclusions have been reached: 

• The scaling limit of g(v), q — > and v ~ q, has the common form for n = 0: 



g(v) ~ q X g(v/q), with A = < 



3, d > 8 and in mean field, 

2-e/2, 6<d<8, 

l + 7 //3 = 2 + e/2 + 0(e 2 ), d < 6; 

(50) 

e = 6 — d. The scaling function has proved to be very robust, as it is 

g{x) ~ x 2 {3 + x) (51) 

for all three cases. 

This result is not unexpected for d > 8, as it must be true for the full g(v) (even for any 
finite n) built up from the exact vertices. This can be proven by applying a simple scaling 
analysis around the Gaussian fixed point. Using (jHJ), (|Cip and (|Q2j) . the exact g(y) can be 
separated into two parts: 

g(y) = (w 2 V 2 + ^U 2 V S + ...)+ ^ I d > 8, U = 0. 

The first term represents the analytical part, the ellipsis dots stand for its subdominant 
contributions with 3 < k < | — 1, whereas the second one is nonanalytical with all of its 
subterms being the same order for q — > and v ~ q: 

£ *. G (k+i) v kc* £ c k <f-t(v/q) k = f-*h(v/q). (52) 
fc>f-i ' fe>|-i 

The analytical part obviously dominates for ci > 8, and the mean field like truncated form 
in Eq. (I24p is valid for the exact g(v) (and even for a generic n too). 



24 



In the transitional regime, 6 < d < 8, the analytical part disappears, and g(y) takes the 
form as in (|52j) . There is now no reason to think that the exact scaling function h has the 
mean field functional form re 2 (3 + x), since all the higher order vertices equally contribute. 
Nevertheless, the first order calculation provides — somewhat surprisingly — Eqs. (|50p and 
(I5ip . suggesting that a cancellation of terms with vertices other than the cubic and quartic 
ones must occur in one-loop order. A similar remark is valid for the regime below the upper 
critical dimension too. It must be realized, however, that the truncated form (|24p or (|51|) 
holds only for the spin glass limit n = 0, as opposed to the d > 8 regime, and an infinite 
number of terms contribute for n finite. 

• The mean field scenario is very robust for decreasing d, see Eqs. (f50|) and (f5Tj) . leading to 
an unstable RS phase for n = even below 6 dimensions. By computing the boundary of 
stability, 2mi = 0, in finite n (a kind of AT line in zero magnetic field), even a shrinking RS 
phase has been found with respect to the mean field picture (Fig. [1]); at least in this first 
order perturbative approximation. 




field, d = oo 



wiq 



FIG. 1: Phase boundary of the RS-RSB transition (Almeida-Thouless line) in the q — n plane. This one- 
loop result shows a shrinking RS phase with decreasing d. (a) and (b) represent (|46p . while the nonanalytic 
feature of the AT line in the transitional regime (6 < d < 8) , Eq. (|47|) , is displayed in (c) . The stable RS 
domain for small but finite n disappears for d < 6 and n<n c ~e = 6 — d, as shown in (d). 

• The mean field formula Eq. (|23p . which connects the replicon mass (2mi) and the replicon 
3-point vertex (u>2) for n = 0, remains valid below 6 dimensions, at least in the leading e 
order. This follows from the robustness of g, as shown in (|50p and (|5ip . From the independent 



calculation of the one-loop 3-point vertices 
u>2 (keeping n finite here): 

4 — n n 
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we can quote the relatively simple result for 



(n-2) 2 ^ v ^' J (2vr) 



(p 2 + n) 2 (p 2 + 2)2 (p 2 + n )3(p2 + 2) 



xe + 0(e 2 ). 
(53) 



Eq. (|49p follows immediately from (|23p in the spin glass limit, providing us a good check of 
the theory. 

As it can be observed from the above discussion, stability of the RS spin glass state in zero 
magnetic field, and thus the possibility of the droplet scenario itself, can not occur in those 
high dimensions where direct perturbation expansion (loop expansion) or the perturbative 
renormalization group is applicable (d < 6). There is even an obvious tendency that RSB 
is escalating with decreasing space dimension. Although common objections against this 
statement might argue that a droplet-like RS state is intrinsically nonperturbative, it must 
be stressed that, at least for d > 8, exact considerations are in accordance with the (first 
order) perturbative result. 

The possibility for the emergence of a marginally stable droplet phase — notwithstanding 
the ever increasing RSB domain — in low enough dimensions will be put forward in this 
remark. The n-dependent scaling form of g(v) for the droplet picture, Eq. (i32|) . proposed 
in Section IIVI clearly shows that g for the droplet model should be a nontrivial generalized 
homogeneous function of v/q and n near T c . Just below 6 dimensions, g of (|48p is not such, 
as it can be decomposed into the form g = g an + g nonan , where g an is analytical (<? n onan is 
nonanalytical) for n -> 0, respectively, and # n onan < San for any v/q. (<? n onan ~ n 1_e/2 
when 6 < d < 8.) As a matter of fact, n in g plays the role of a small mass, the "replicon" 
mass, as contrasted to the large mass, the "longitudinal" one. This small mass is, however, 
innocuous around d = 6, since it does not cause infrared problems. (See as an example Eq. 
f)53[) : an infrared divergence would be expected at d = 4, if the collapse of the perturbative 
approach were neglected.) We can postulate a dimension cZr, which is only naively believed 
to be 4, where </ n onan stops being negligible, and g may become a nontrivial generalized 
homogeneous function of v/q and n. Although a proof of this scenario and an estimate of 
c?r seem to be a rather hard task — due to the nonperturbative nature of the problem in 
these low dimensions — , these considerations could help the understanding how a marginally 
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stable replica symmetric phase may emerge in low enough dimensions. 5 



32l | : In this paper the authors use the 



Our final comment concerns the work of Bray and Moore 
same "two-packet" RSB for temporarily breaking the replica symmetry, and eventually restoring 
the RS state by sending p (which they call m) to infinity. What is more, this RS state is marginally 
stable! It must be stressed, however, that the two approaches differ substantially: p < n is always 
assumed here, and the "two-packet" RSB calculation used in this paper is only a tool to compute 



32] is 



the replicon generator of the RS state. The marginally stable RS phase found in Ref. 
obviously in conflict with our results, since the replica symmetric spin glass state has proved to be 
unstable in the regime 6 < d < 8, assuming only the invariance of the Legendre-transformed free 
energy, Eq. (|7|), as a prerequisite. 



APPENDIX A: BASIC DEFINITIONS 



In this appendix the definitions of the lower and upper case vertices are displayed in a tabular 
form for the generic RS state. Although not all of them are used in this paper, completeness was 
favoured up to the quartic order. Due to the RS symmetry, the number of different vertices is 
limited to 3 for the masses (k = 2), 8 for the cubic (k = 3) and 23 for the quartic (k = 4) ones. 
The first column provides an obvious graphical representation, while in the second one the lower 
case bare couplings are enumerated together with their corresponding invariants; see Eqs. (|4|) and 
([5]). The alternative set of vertices with the upper case notation is best defined by the derivatives 
of the exact Legendre-transformed free energy F, Eq. (I10p . and they are listed in the last column. 
The corresponding bare couplings enter the Lagrangian in the following way (see footnote H]): 

C = \T. E {p 2 SKl lS + M aPaS )tfcf>l S p +£\ (AI) 

P (a/3),( 7 <5) 

with the interaction part 

(ct/3) k=3 fc - iV2 V^ (ai/3!),...,(a fc /3 fc ) 

[H = 2h follows trivially by comparing ([5]) and (|A2|) .] 

As already declared in the main text, the common notations are used in what follows for the 
quadratic {—fhj and —Mj), cubic (Wj and Wj) and quartic (Uj and Uj) vertices (instead of g)j 
and Gf\ k = 2,3 and 4). 



dn, the lower critical dimension of RSB, should not be missed with di, the lower critical dimension of the spin 
glass state. Of course, the transition from RSB to the droplet scenario can occur only if d; < dn- 



The masses, k = 2: 



j 




3 




1 




a/3 


-Mi = -M^^ 


2 




a/37 


-M 2 = -M Q7)/ 3 7 


3 


• • 


a/375 


-M 3 = -M a/3)7<5 



Cubic vertices, k = 3: 



J 








1 


A 


^E«c 

a/37 




2 




a/3 


^2 = ^0/3,0/3,0/3 


3 




^E<*« 

a/37 


W3 = W Q/ 3 iQ/ 3 ij( 3 7 


4 


<z> • • 


^ E <*gM 

a/37<5 


W 4 = W Q/ 3,a/3,75 


5 




^ e <c< 

a/37<5 


W5 = W Q/ 3 iQ7j/ 35 


6 


A 


^ E <«f 

a/375 


^6 = W a p ja ~f >a s 


7 




^ E C«<* 




8 


• — • 
• — • 


E «c 

a/375/11; 


W 8 = W a p >y s,nv 



Quartic vertices, k = 4: 
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1 

J 




=(*) /■(*) 


G (fc) 


1 


I 1 

i 1 


mi V 6 al3 d, 13 ' 1 cP s (b 5a 

Q/375 


C^l = ^a/3,/37,75,5a 


2 






^2 = ^a/3,a/3,a/3,a/3 


3 


<=xz> 


a/37 


^3 = ^a7,a7,/37,/37 


4 


<z> <z> 


a/375 


[^4 = ^a/3,a/3,75,75 


5 


A 


a/37 


^5 = ^a/3,a/3,a7,/37 


6 


>~ 


E <c*£>£ 

a/37<5 


= ^a/3, 07,^7, 75 


7 


> — 


-7 E 

a/375/1 


^7 = Uafoayfl^Jn 


8 


«o ■ 


E 

a/37 


Us = ^a/3,a/3,a/3,/37 


9 




-9 E <*S?*SM 
a/375 


= ^a/3,a/3,a/3,75 


10 


<=< 

>• 


«io E 

a/375 


C^IO = ^a/3,a/3,/37,/35 


11 


o ■ ■ 


V 6 al3 6 aP dJ 3 ^ dp 5 

a/375 


^11 = f^a/3,a/3,/37,75 


12 




«« E 

a/375/i 


C^12 = ^a/3,a/3,/37,5/t 


13 


...^ 


«is E <«c« 
a/375 


^13 = ^a/3,a/3,a7,/35 


14 


>• 


-14 E 

a/375/i 


f 14 = ^a/3,a/3,75,7/i 


15 




«« E 

a/375/11/ 




16 




-16 E <«^3* 
a/375/i 




17 




Sir E tf^tfsC 

a/3 , ySiM> 


Un = U a /3 : py tl s,iJ,v 


18 


~< 
>• 


-18 E «1^3< 
a/375/i 


Uis = U a p^ n s,^n 


19 


x 


-19 E «c< 


U\g = f^ta,/i/3,/t7,/i5 


20 


- <• 

>• 


-20 E 

afl'ySfj.i/ 


^20 = £^0/3,75,7/4,71/ 


21 


> <. 


"^1 Z-/ ^Pl^P2^P3^P4 

a/375^1/ 




22 




-22 E 

a/375/i^p 


U22 = Uay^Swp 


23 




-23 E «cc 


U23 = U a p n s,tiv,pu 



29 



APPENDIX B: THE LINEAR RELATIONS BETWEEN THE TWO SETS OF VERTICES 



Decomposing the restricted sums in Eq. (1A2I) and comparison with (JH) and © provide the 
linear relations between the upper and lower case bare couplings. (Although this is straightforward 
combinatorics, it becomes quite time consuming for highly disconnected quartic vertices.) In what 
follows, these formulae are displayed by omitting the bars in the notations; in this form they can 
be considered as the definition of the exact lower case vertices. The results for the masses and 
cubic vertices are taken over from 



mi = -(Mi-2M 2 + M 3 ) 
m 2 = M 2 - M 3 
m 3 = jM 3 ; 



431 ] . and they are shown here for easing to refer them: 
Wl = W 1 - 3W 5 + 3W 7 - W 8 

w 2 

W3 = :}]} - :;1T| - nil - 3TT',; + 1511V - 01! s 

w i 

3W 5 - 6W 7 + 3Wt 
W 6 - 3W 7 + 2W 8 



^W 2 - 3W 3 + ^W 4 + 3W 5 + 2W 6 - 6W 7 + 2W 8 



(Bl) 



w 5 

Wq 
W 7 



The set of equations for the quartic case: 



\w 7 - ^W 8 
\w 8 . 



Ul = 3U! - 12U W + 12U 17 + 6U 2l - 12U 22 + 3U 23 

u 2 = 3U X + -U 2 - 3U 3 + |z/ 4 - 4C/ 8 + 2U 9 + 12U 10 + 12U n - 24U 12 + 6U 13 - 12U U + 12U 15 
-24U W + 48U 17 - 24U 18 - 6U 19 + 24U 20 + 30U 21 - 72U 22 + 18U 23 

u 3 = -6Ui + 3U 3 - 3t7 4 - 6Uio - 12U U + 12Ui 2 + 18U M - 121% + 36U W - 48Ui 7 + l2Ui 8 
+3U 19 - 12U 20 - 39U 21 + 72U 22 - 18U 23 

Ui = hj A - 3U U + ~U 15 + 3U 2l - 3U 22 + hj 23 

u b = QU 5 - 24U 6 + 12U 7 - 12U U + 12U 12 - 6U 13 + 6U U - 6U 15 + 36U m - 84U 17 + 48U 18 

-24U 20 - 36U 21 + 96U 22 - 24U 23 
u & = 12U 6 - 12U 7 - 24U 16 + 72U 17 - 24U 18 + 12U 20 + 24U 21 - 84U 22 + 24U 23 
u 7 = 2U 7 - 6U 17 + 6U 22 - 2U 23 

u 8 = 4U 8 - 4U 9 - 12U W - 12Un + 48U 12 - 12U 13 + 12l^i 4 - 24U 15 + 24U 16 - 96U 17 + 48f7i 8 
+8U l9 - 56f/ 20 - 48U 21 + 168C/ 22 - 48U 23 
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u 9 = U 9 - 6U l2 + 3U 15 + 6U l7 + AU 20 ~ 12C/ 22 + 4C/ 23 

-"10 = QU W - 12U 12 - 6U 14 + I2C/15 - 12C/i6 + 48U 17 - 12E7i 8 - 6U 19 + 2AU 20 + 30*7 2 i - 108C/ 22 
+36t7 23 

till = 12E7ii - 12t/i2 - 12U U + I2C/15 - 24C/ 16 + 48^i 7 - 12i7 18 + 12U 20 + 36U 21 - 84U 22 + 24U 23 
«12 = QU 12 - 6C/15 - 12?7i7 - 6C/ 20 + 30f7 2 2 - 12U 23 

ni 3 = -12C/ 12 + 6C/13 + 6C/15 + 36C/ 17 - 24U 18 + 24C/ 20 + 12U 21 - 72U 22 + 24C/ 23 

u u = 3U U - 3U 15 - QU 21 + 9U 22 - 3U 23 

3 9 
U15 = JU15- 6U 22 + -U 23 

u 16 = 12U 16 - 24U 17 - 12U 21 + 36C/ 22 - 12U 23 
un = 6U 17 - 12U 22 + 6U 23 

uis = -24£/ 17 + 12E/i8 - 12Z7 20 - 12C/ 2 i + 60f/ 22 - 24[/ 23 

uig = CTig - 4C/ 20 - 3U 21 + 12U 22 - 6U 23 

u 2 o = 2U 2 o — 6U 22 + 4l7 23 

U21 = 3C/21 — 6f/ 22 + 3C/ 2 3 

3 3 
^22 = 2 22 ~~ 2^ 23 

1 

1*23 = 7^23- 

lb 

APPENDIX C: THE COEFFICIENTS, k = 2, 3 AND 4, IN TERMS OF THE 

VERTICES 

Eqs. Q and ([TO]) define the GW's, which build up 5(f). [See also the paragraph below ([TO]) .] 
For fe = 2, it was worked out in Section HTJ see p4|) and (|15j) : 

—G^ = [2m 1 + n m 2 ] + p(n — p) 4m3. 

We only quote the results for fc = 3 and = 4: 

G (3) = [2u> 2 + nw 3 + n 2 w e ] +p(n-p) [(4iu 4 + 2u> 5 - 2u> 6 ) + n2w 7 ] + p 2 (n - p) 2 8w 8 (CI) 
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and 

G^' = [2u2 + n(u 3 + Us) + n 2 ui + n 3 uig] + p(n - p) [(2ui + 4u 4 + 4u 9 - 2ui + 2u\i + 2ui 3 ) 
+n (2ui2 + 2ui4 + ni6 + «is - 3uig) + n 2 (2u 20 + «2i)] (C2) 
+p 2 (n -pf [(8U15 + 4ui 7 - 4u 20 + 24u 22 + 384u 23 ) +n4u 22 ] +p 3 (n - pf 16u 23 . 

From these examples, we can observe the following properties valid also for generic k: 

is a double polinomial of p[n — p) and n. We must, however, notice that the vertices 
themselves depend on n. When constructing g(v) by its definition ([8]), we are particularly 
interested in the limit p(n — p) — * while n fixed, and then taking the spin glass limit n — > 0. 

• This limit results the fully replicon vertices 2mx, 2u; 2 , 2it 2 , . . . etc. (These are represented 
by simple two-node graphs connected by k edges. See the graphical representations in the 
second column of the tables in Appendix El) This justifies, at least in the spin glass limit, 
the name replicon generator for g(v). 

APPENDIX D: THE LEGENDRE-TRANSFORM METHOD TO DERIVE EQ. <[38jl 

The Legendre-transformed free energy F(v a g) for the model of Eqs. (|Aip and (|A2p . or equiva- 
lently (HJ) and ([5]), is defined by the usual rules 

-F = \nZ-NJ2H a/3 v afi and = Nv aP , (Dl) 

(a/3) aP 

where the partition function Z = J T>cf)e~ l ~ acquires its dependence on the H a p y s by adding the 
source term —N2 Y^( a p) Ha/3 4>p=o to the Lagrangian. By the second equation of (|D1[) y/N v a p = 
(0p=o) f°h° ws (the average is taken by the measure e _-C ), and by the following redefinition of the 
fields we can make the 1-point function disappear: 

<pf = cl>f-VNv aP d^ =0 , resulting (<ptf) = 0. (D2) 

The reader is now reminded that the bare couplings of the original sourceless — i.e. H a p = - 
system are defined likewise, (<^p£ ) = [see the remark in the paragraph below ((3J)]. Hence we 
can see that v a p = if H a p = 0, and it follows from the inverse expression ^ = NH a p that 
stationarity of F occurs for v a p = 0. 
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For computing F, (f> is replaced by ip, according to Eq. (|D2h . in the Lagrangian [ (|Alj) and (|A2|) . 
supplemented by the source] to give from (|D1|) 

+ H^v aP -- M a p jlS v aP v lS + X) y Yl GaiP 1 ,..., ak fo v *iPi--- v *A- ( D3 ) 

(a/3) (a/3), ( 7 <5) fc=3 ' (ai/3 1 ) ) ...,(a fcJ 9 fc ) 

The new Lagrangian £ of the shifted system gives In Z, its coupling parameters are listed below 
up to cubic order: 

Ha/3 = (H + H a p) - ^ M a p a s VjS + ^ ^ W a /3,jS,pu V-ygVpu 
(7(5) (7c?), (pV) 



+ ..., 



6 

( 7 <5),(p^),(pw) 



1 - (D4) 

M a p n S = M a p n s — Waf3,~r6,fiv V^ u — - U a p n S,pv,pu) V^Vpuj + . . . , 

(pu) (pu),(pu) 

^^a(3,^S,pu — W^af3,y8,pv ~i~ ^ ^ £^a/3,7<5,p^,pu; ^ pu; ""(-••• 

(pw) 

Since -£f aj g enters C solely via H and the condition in (|D2p guarantees the lack of 1-point insertions, 
F depends on the v a ps alone. 9 

As it is well known from the literature (see (44! \^\ , for instance, and references therein), F is 
the generator function for the exact vertices, and since 

F( v af3) = F(q + v ap) + [(/-dependent term], 
[cf. with CED], we have by Eq. CE0|): 

-^=llnZ(^ = 0) 

1 1 

- 2 2^ M ^.7<5 ^/3^ + fci ZZ G ai(3 lt ...,a h p k <Wr ■ ■ ■ (D5) 

(a/9),(7<5) ^=3 ' (ai/3i),...,(Q fc /3 fc ) 

Proceeding further along the same lines with Section IIII] the one- variable function F(v) is defined 
by taking v a g = v for a n (3 = 0, otherwise u Q/ g is 0. From the two alternative forms for F in Eqs. 
(|D3P and (|D5p . and using the generic definitions (JSj) and ©, our basic result for g(v), displayed in 
Eq. (J3ll), follows. 



9 Although the derivation of F presented here was for the RS spin glass field theory, the method is entirely general: 
instead of (a/3) , a generic index set for the field components could have been used, and the symmetry of the model 
does not effectively matter. 
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APPENDIX E: BLOCK DIAGONALIZING A MATRIX WITH THE "TWO-PACKET" 

REPLICA SYMMETRY BREAKING 



In this Appendix the greek-roman notation of Bray and Moore 32] is used to ease the classi- 
fication of matrix elements when replica symmetry is broken by dividing the n replicas into two 
groups: the p replicas in the first group are denoted by a, /3, . . . , while the other n — p by a, b, 
. . . . For the 15 different matrix elements the following concise notations will be used: 

• First class (diagonal elements): 

(U = M a g a g 

W — H/T 



M{> = M, 



M[ {2) =M ab , ab 



Second class (one common replica): 



M 2 (1) = M a7ife7 M 2 (2) = M a7i/37 M 2 (3) = M a7jfe7 

M^ (1) = M ac ^ c M'l 2) = M acM Mf ] = M ac , bc 



• Third class (all replicas are different): 

m M 3 (2) = M Q/3)7 5 M 3 (3) = M af3 ^ d 

M^> = M aPAb M 3 W = M abnd 

M^> = M abyCd = M ab)Cd 

The procedure of Ref. [48] can be applied to find the block diagonalized matrix M by a similarity 
transformation (see also [43]]). The basis vectors 10 of the invariant subspaces, matrix blocks and 
multiplicities are listed below: 

1. The 3-dimensional longitudinal subspace: 

*(L0) = 1 <1) = 1 ^(Ll') = 1 



Any elements of the basis vectors not indicated explicitly are zero. 
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M ( lo),(lo) = + (n - V ~ + (P~ + (p - l)(n - p - 1)M 3 (4) 

M (L0 ),(li) = ^(p - 1) [(P - 2)M 3 (3) + 2M 2 (3) ] 
M (LO ),(LiO = \{n-p-l)[{n-p- 2)i\4 (3) + 2M^ (3) ] 



%i),(LO) = {n-p) [2M 2 (3) + (p - 2)M 3 (3) ] 

%i),(Li) = M[ 2) + 2(p - 2)M 2 (2) + i(p - 2)(p - 3)M 3 ( 



M(Li),(Li') = ~(n - P)(n - P - 1)M 3 (1) 



M mm = P [2M^ (3) + (n - p - 2)M^ (3) ] 

M( L 1'),(L1) = ^1>(P " 1)^3 (1) 

M (L iO,(LiO = M[ {2) + 2(n - p - 2)M^ (2) + i(n - p - 2)(n - p - 3)M^ (2) 



Multiplicity: 3 



2. The 2-dimensional anomalous (A) subspaces: 



Co) 



1 if a /i 

— (p — 1) if a = u 



V) 



1 if a, /? n 

— 7j(p — 2) if a = /x or (3 = n 



Vo) 



1 if a 7^ m 

— (n—p — 1) if a = m 



1 if a, 6 m 

- i (n — p — 2) if a = m or b = m 



%o),(M)) = M[ 1] + (n-p- 1)M 2 (1) - M' 2 {1) - (n - p - 1)M 3 (4) 



p-block: 



%o)^i) = ^(p-2) (M 2 (3) -Mf) 



r(3)> 



%i )l(At o)=2(n-p) (M 2 (3) -Mf) 

= M[ 2) + (p - 4)M 2 (2) - (p - 3)M 3 (2) 



r(3)> 



m-block: 



M, 



(mO),(ml) 



M, 



(ml),(mO) 



M(mO),(mO) = + (p - 1)M^ 1} - M 2 « - (p - 1)M. 



(1) 



r( 4 ) 



-(n-p-2) (M^-M^) 



r/(3)> 



2p(M 2 (3) -M^) 



/(3)> 



M (ml) , (ml) = M{ (2) + (n - p - 4)M^ - (n - p - 3)M£ 



/(2) 



r/(2) 



Multiplicity: 2 [(p - 1) + (n - p - 1)] = 2(n - 2) 



The 1-dimensional RO replicon subspaces: 



1 if a / p and a ^ m 

— (p — 1) if a = p and a ^ m 

— (n — p— 1) if a 7^ p and a = m 
{jp — 1) (n — p — 1) if a = p and a = m 



Multiplicity: (p — 1) (n — p — 1) = p(n — p) — (n — 1) 
The 1-dimensional Rl replicon subspaces: 
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• (fj,u) eigenvector: 
1 



V) 



i(p-3) 



a, (3 ^= fx, v 

if one of a and /3 coincides with one of and z/ 



5 (p — 2) (p — 3) if a = n and /? = u, or a = v and (3 = /J, 



M ( ^),(^)=M 1 (2) -2M 2 (2) +M 3 (2) 



Multiplicity: -p(p — 3) 



(mn) eigenvector: 



-i(n-p-3) 



if a, b / m, n 

if one of a and b coincides with one of m and n 



i (n — p — 2) (n — p — 3) if a = m and 6 = n, or a = n and b = m 



M 



(mn),{mn) 



M[ {2) -2M^' Z> + 



'( 2 ) , ^'(2) 



Multiplicity: - (n — p) (n — p — 3) 
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